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ABSTRACT 

Equations  governing  the  motion  of  a  specific  class  of  singularities  of  the  Euler  equation  in 
the  extended  complex  spatial  domain  are  derived.  Under  some  assumptions,  it  is  shown  how 
this  motion  is  dictated  by  the  smooth  part  of  the  complex  velocity  at  a  singular  point  in  the 
unphysical  domain.  These  results  are  used  to  relate  the  motion  of  complex  singularities  to 
the  stability  of  steady  solutions  of  the  Euler  equation.  A  sufficient  condition  for  instability 
is  conjectured.  Several  examples  are  presented  to  demonstrate  the  efficacy  of  this  sufficient 
condition  which  include  the  class  of  elliptical  flows  and  the  Kelvin-Stuart  Cat’s  Eye. 
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1.  INTRODUCTION 


Singularities  of  the  Euler  equation  for  inviscid  incompressible  flow  have  been  studied  ex¬ 
tensively  for  many  years.  An  important  question  that  remains  quite  controversial  is  whether, 
in  three  dimensions,  an  Euler  flow  can  develop  a  singularity  in  finite  time  when  the  initial 
conditions  are  smooth  (see  Pumir  &:  Siggia1,  Majda2,  Kerr3,  Brachet  et  al4  and  the  references 
therein).  While  this  paper  does  not  address  this  controversy  directly,  except  within  a  narrow 
class  of  singularities  under  restrictive  assumptions,  we  expect  that  the  approach  presented 
here  may  help  in  a  fuller  resolution  of  this  problem.  All  previous  singularity  studies  have 
been  limited  to  examining  solutions  to  the  Euler  equation  in  the  real  (physical)  spatial  do¬ 
main,  though  in  some  cases  the  singularity  presence  in  the  complex  unphysical  plane  was 
deduced  from  monitoring  the  Fourier  coefficients. 

For  simpler  Euler  flows  with  sharp  discontinuities  like  the  Kelvin- Helmholtz  problem, 
previous  work  by  Moore5,6,  Caflisch  &.  Orletma7  and  Cowley  et  al8  have  shown  the  advan¬ 
tage  of  studying  the  equations  in  the  complex  unphysical  plane.  In  that  case,  the  singularity 
formation  in  the  physical  domain  can  be  understood  in  terms  of  two  processes:  (a)  instanta¬ 
neous  singularity  formation  at  some  points  in  the  complex  unphysical  domain  at  the  initial 
time,  and  (b)  the  motion  of  such  complex  singularities  towards  the  physical  domain.  The 
question  of  finite  time  singularities  in  the  real  domain  then  becomes  a  question  of  whether 
such  singularities  approaching  the  real  domain  actually  impinge  it  in  finite  time.  Moore5,6 
and  Caflisch  &  Orlenna7  have  shown  that  for  specific  symmetric  initial  conditions,  the  sin¬ 
gularities  form  on  the  imaginary  axis  after  a  finite  time,  which  then  move  towards  the  origin 
leading  to  a  singularity  formation  on  the  real  axis.  However,  recent  work  by  Cowley  et  al8 
suggests  that  the  Moore  singularity5,6  actually  forms  at  the  initial  instant  of  time  at  points 
off  the  imaginary  axis.  Later  on,  some  of  these  singularities  merge  in  pairs  on  the  imaginary 
axis  before  traveling  towards  the  origin  (which  is  part  of  the  real  domain).  Indeed,  it  appears 
that  in  a  wide  class  of  differential  equations  including  the  familiar  one-dimensional  inviscid 
Burger’s  equation,  singularities  that  show  up  later  in  the  real  domain  are  either  present  at 
some  complex  location  initially  or  are  formed  at  the  initial  instant  of  time.  Such  singulari¬ 
ties  typically  form  at  points  (including  infinity)  in  the  complex  spatial  plane  where  a  regular 
perturbation  series  in  time  (for  early  times)  becomes  disordered.  In  some  cases,  it  is  found8,9 
that  singularities  are  born  at  infinity  in  the  sense  that  the  singularity  location  in  the  finite 
spatial  plane  recedes  to  infinity  as  t  — >  0+.  The  study  of  complex  singularities  appears  to 
have  been  a  fruitful  exercise  that  transcends  the  question  of  finite  time  singularity.  Indeed, 
recent  work  on  an  interfacial  evolution  problem  by  Tanveer10  suggests  that  this  may  be  an 
important  step  in  understanding  the  effect  of  small  regularization  on  an  otherwise  ill-posed 
time-evolving  flow  in  the  physical  domain. 
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Despite  successes  in  each  of  the  cases  cited  above,  and  perhaps  in  other  cases  unknown 
to  the  authors,  the  problems  were  reducible  to  one  independent  space-like  variable  besides 
time.  We  are  unaware  of  attempts  to  generalize  these  ideas  for  multiple  spatial  variables 
as  appears  necessary  for  studying  a  general  Euler  flow.  Here,  we  study  a  class  of  moving 
singularities  of  the  Euler  equation  whose  form  appears  to  remain  unchanged  with  time,  with 
the  expectation  that,  for  certain  other  initial  conditions,  the  process  (a)  described  in  the  last 
paragraph  leads  to  the  formation  of  complex  singularities  of  this  form  when  they  are  not 
present  initially. 

In  this  paper,  we  provide  analytical  evidence  to  show  how  each  singularity  of  a  specific 
type  moves  with  the  smooth  part  of  the  complex  velocity  field  evaluated  at  that  point.  Of 
course,  in  a  simpler  vein,  it  has  long  been  known  that  point  vortex  singularities  in  the  physical 
domain  move  with  the  local  velocity  at  the  singularity  location11.  However,  with  certain 
assumptions,  this  result  is  generalizable  to  singularity  motion  in  the  complex  unphysical 
plane  as  well.  This  result  is  used  to  conjecture  sufficient  conditions  for  the  instability  of  a 
steady  solution  of  the  Euler  equation.  We  present  two  conjectures.  Conjecture  1  addresses 
the  absolute  instability  of  any  steady  Euler  flow  with  a  stagnation  point.  Conjecture  2 
addresses  instability  in  the  Langrangian  sense;  when  the  conditions  of  this  conjecture  hold, 
there  exist  disturbances  that  grow  as  we  follow  some  fluid  particle. 

Several  applications  to  the  hydrodynamic  stability  of  inviscid  flows  in  an  infinite  or  peri¬ 
odic  domain  -  which  include  the  class  of  elliptical  flows  and  the  Kelvin-Stuart  Cat’s  Eye  - 
are  made  to  illustrate  the  use  of  Conjecture  1  to  determine  the  absolute  instability  of  these 
flows. 

2.  MATHEMATICAL  EQUATIONS  IN  THE  COMPLEXIFIED  DOMAIN 

Consider  the  Euler  equation  for  the  motion  of  an  invisicid  fluid  where  the  vorticity  uj(x,  t ) 
and  the  velocity  u(x,t)  are  determined  by 

ut  +  u  •  Vu;  —  u  ■  Vu  =  0  (1) 

V  •  u  =  0  (2) 

V  x«  =  w  (3) 

Each  of  the  components  of  the  above  equations  can  be  put  in  the  form 

F(xux2,x3,t)  =  0  (4) 

where  x  =  (x\,  x2,  x3).  In  the  case  of  27T  periodicity  in  each  direction,  if  we  assume  that  each 
coefficient  C]ki  in  the  following  representation  of  F 

F(x i,  x2,  x3)  =  £  Cjkl  e«l»*.+**a+'*3]  (5) 

j,k,l 
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satisfies  the  relation 

\Cjk,\  <  K  e-5(|j|+|fc|+|/|)  (6) 

for  some  constant  K  and  8  which  are  independent  of  j,  k  and  /,  then  it  follows  that  (4)  will 
hold  for  complex  ®i,  x 2  and  £3  when  |/m  £i|,  \Im  £2!  and  \Im  £3!  are  each  small.  Since  the 
right  hand  side  of  (4)  is  zero,  which  is  an  analytic  function  of  x\,  £2  and  £3  in  the  complex 
domain,  it  follows  that  (4)  holds  for  all  complex  x.  The  condition  (6)  will  certainly  hold  for 
the  Fourier  series  coefficient  of  each  component  of  equation  (1),  (2)  and  (3)  if  each  component 
of  u  and  u  and  its  derivatives  satisfy  such  a  condition.  We  restrict  our  attention  to  cases 
when  this  is  true  at  some  early  time,  which  we  assume  is  possible  for  a  certain  subclass  of 
initial  conditions.  Then  we  arrive  at  the  conclusion  that  each  of  equations  (1),  (2)  and  (3) 
hold  for  complex  £  (i.e.  each  component  is  complex  ).  However,  since  the  right  hand  side 
of  each  component  of  (1),  (2)  and  (3)  is  zero,  which  is  analytic  for  all  complex  x  and  t,  it 
follows  that  (1),  (2)  and  (3)  will  hold  for  all  x  and  t.  The  singularities  of  u  and  u>  are  those 

that  are  consistent  with  the  solutions  of  (1),  (2)  and  (3)  in  the  complex  x  domain.  The 

same  arguments  hold  even  when  the  periods  in  x\  ,  x2  and  £3  are  not  the  same  or  when  the 
geometry  is  infinite  in  some  directions  and  periodic  in  others. 

We  examine  a  certain  kind  of  initial  conditions  for  which  the  velocity  and  vorticity  are 
each  real  and  analytic  everywhere  for  real  £,  but  have  complex  singularities  at  a  singularity 
surface  determined  by 

d{x,0)  =  0  (7) 

where  d  is  an  analytic  function  of  each  component  of  x  that  is  real  valued  and  greater  than  0 
for  real  x.  Here  d  is  also  a  function  of  t  and  its  time  evolution  equation  will  be  derived  later. 
We  also  require  the  non- degeneracy  condition  that  Vd  /  0  where  d  =  0*.  We  consider 
initial  conditions  that  can  be  decomposed  into 

u(x,  0)  =  u,(x,  0)  +  f(d(x,  0))  q{x ,  0)  (8) 

w(£,  0)  =  us(x,  0)  +  f'(d(x,  0))  p(£,  0)  (9) 

where 

f(d)  =  d°  (10) 

with  0  <  a  <  1.  Thus  the  initial  vorticity  blows  up  at  thi  complex  singular  points,  though 
the  velocity  is  finite.  In  (8)  and  (9),  we  assume  that  f.ich  of  u3(x,  0),  ujs(x,  0),  q(x,  0)  and 
p(x,  0)  are  analytic  functions  of  x  ana  that  they  are  real  valued  for  real  x.  Thus  the  initial 

*We  have  shown  later  that  this  condition  is  satisfied  at  later  times  once  satisfied  initially. 
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conditions  (8)  and  (9)  correspond  to  smooth  initial  conditions  in  the  real  domain.  We  also 
require  that 


V  •  us(:r,  0)  =  0 

("I 

V  x  us(x,  0)  =  u3{x,  0) 

(12) 

and 

V  x  [/  g(x,0)]  =  f'p{x,0) 

(13) 

Now,  for  t  >  0,  we  seek  solutions  to  (1)  -  (3)  of  the  form 

u(x,  t)  =  us(x,  t)  +  f{d(x,  t ))  q(x,t) 

(14) 

u(x,  t )  =  LJs(x,  t)  +  f\d{x,  t ))  p(x,  t) 

(15) 

From  (2)  and  (3),  we  get  the  following  relations: 

T  f 

V  •«,  +  /'  ~S/-q  +  q-\/d  =  0 

(16) 

r/- 

V  x  +  /'  x  q  +  Vd  x  q  =  us  +  fp 

(17) 

If  we  now  require  that 

V  •  u3  =  0 

(18) 

and 

V  x  u  s  =  u3 

(19) 

wS|  +  u,  •  VcDs  -  u>,  •  Vus  =  0 


then  the  relations  (16)  and  (17)  imply  that 


Lv-q  +  q-Vd=  0  (21) 

X  q  +  Vdxq=p  (22) 

In  addition,  the  substitution  of  (14)  and  (15)  into  (1),  with  the  use  of  relations  (20)  -  (22), 
yields 

f'W,  +  a, •  Vrf]  +  f  p,  +  a,  ■  Vp-p^-v ■  q  +  /(,-•  V)  p  +  Lq-  V)  a, 

■  v  d)  -  (p  •  v)  s,  -  L  (a. .  v)  ?/(W)  •  (v  x  ,-)  -  if-  v?  =  o  (23) 
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Now  notice  that  if  we  choose  to  evolve  d(x,t)  according  to 


dt  +  us-Vd  =  0  (24) 

then  it 

Pt  +  us  •  Vp  -  p^p-V  '  9  +  /(<7  •  V)  p  +  —(</*•  V)  u>s  =  q(us  •  V  d) 

+  (p  •  V)  +  j,  (us  ■  V)  9  +  <f/(Vd)  •  (V  x  (f)  +  /p  •  V<f  (25) 

The  existence  of  p  satisfying  (21),  (22)  and  (25)  for  small  enough  time  follows  from  the 
existence  of  u  and  wls,  each  satisfying  (1)  -  (3),  and  of  d  satisfying  (24),  since  u  —  us  =  f'{d)p. 
We  also  notice  that  the  coefficients  in  (21), (22)  and  (25)  are  continuous  even  at  d  -  0. 
Therefore,  it  appears  reasonable  to  assume  that  for  small  enough  time,  with  appropriately 
regular  initial  data,  such  a  solution  is  C1  in  Xj,  x-i  and  X3  at  d  —  0  in  the  extended  complex 

domain.  Since,  (24)  implies  that  d  >  0  in  the  real  domain,  once  it  is  so  initially,  it  follows 

that  such  a  solution  p  will  be  real  and  analytic  in  the  real  domain  for  sufficiently  small  t. 
since  it  is  so  initially.  Notice  that  if  we  were  to  choose  dt  +  us  •  Vd  ^  0  at  d  =  0,  then  the 
resulting  evolution  equation  for  p  would  contain  singular  coefficients  in  x  at  the  singularity 
surface  d  —  0,  which  would  preclude  a  C1  solution  for  p  at  d  =  0.  Indeed  for  a  general  choice, 
p  is  not  even  expected  to  be  finite  at  d  =  0. 

Despite  the  choice  (24)  for  the  evolution  of  d  -  which  results  in  (25)  having  C1  coefficients 
when  d  =  0  -  we  cannot  rule  out  the  possibility  that  the  system  of  nonlinear  evolution 
equations  for  cDs,  p  and  d  develops  spontaneous  singularities  after  a  finite  time.  Hence, 
global  existence  in  time  cannot  be  expected  in  general.  However,  we  assume  that  there 
exists  a  nontrivial  subclass  of  initial  conditions  so  that  no  spontaneous  singularity  develops 
in  the  complex  domain  for  p,  q,  u3  or  us,  other  than  a  weak  singularity  of  p  and  q  at  d  =  0 
(in  the  sense  that  p  is  Cx  but  not  analytic).  In  that  case,  the  continuity  of  p  at  d  =  0  implies 
that  at  later  times,  the  singularity  of  u  is  still  determined  by 

d(x,*)  =  0  (26) 

where  /'  is  singular.  Writing  (24)  in  the  characteristic  form,  it  is  clear  that  (26)  will  be 
satisfied  for  points  x0{t)  determined  by 

=  (27) 

with  the  initial  condition  xo(0)  satisfying  the  constraint  that 

d(xo(0),0)  =  0  (28) 
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Thus,  each  point  on  the  singulaiity  surface  moves  with  the  speed  us  evaluated  at  that  point. 
However,  us  is  the  smooth  part  of  the  vector  field  u  at  d  =  0  in  the  complex  x  plane,  which 
means  the  location  of  complex  singularities  of  the  type  considered  here  can  be  determined 
without  having  to  calculate  p(xj)  from  (25).  Furthermore,  with  the  assumption  on  the 
smoothness  of  us  and  continuity  of  p  for  all  time,  any  singularity  of  this  kind  can  only 
approach  the  physical  domain  without  actually  impinging  it  in  finite  time.  This  can  be 
deduced  from  (24)  by  putting  it  in  the  characteristic  form  and  noting  that  us  is  real  in  the 
real  domain  and  d(x,  0)  >  0  in  the  real  physical  domain.  Furthermore,  if  a  flow  has  a  real 
point  (say  the  origin,  without  loss  of  generality)  where  us  =  0  with  a  velocity  gradient  tensor 
l\j  =  independent  of  time,  then  replacing  us  in  (27)  by  T  x  shows  that  a  complex 
singularity  can  only  approach  the  physical  domain  exponentially  in  time.  This  observation 
may  be  relevant  to  the  numerical  calculation  of  Brachet  et  al4  in  the  Taylor-Green  problem, 
where  they  note  an  exponential  approach  of  a  singularity  towards  the  physical  domain. 
However,  since  the  deformation  tensor  is  traceless,  because  of  the  divergence  condition,  it  is 
clear  that  only  a  restricted  set  of  points  on  the  initial  singularity  surface  are  a  subset  of  the 
complex  stable  manifold  of  the  fixed  point5. 

Equation  (27)  is  convenient  since  it  shows  that  one  need  not  solve  for  p(x,t)  to  determine 
the  singularity  location  in  the  complex  plane  even  though  one  needs  to  determine  us(x,t) 
in  the  complex  plane  When  exact  solutions  are  known  for  us ,  which  may  be  steady  or 
unsteady,  this  poses  no  problems;  however,  when  the  solution  is  only  determined  numerically 
in  the  real  domain,  one  cannot  expect  to  determine  us  in  the  complex  plane  through  Fourier 
extrapolation,  since  such  a  procedure  is  ill-posed.  In  such  a  case,  the  solutions  have  to  be 
determined  numerically,  directly  in  the  complex  plane.  This  will  be  the  subject  of  another 
paper.  However,  if  the  singularities  are  close  to  the  physical  domain,  it  is  possible  to  obtain 
complex  singularity  trajectories  directly  from  us  in  the  real  domain.  This  can  be  seen  by 
decomposing  Xo{t)  =  x(t)  +  iy(t),  into  its  real  and  imaginary  parts  and  taking  the  real  and 
imaginary  part  of  (27)  in  the  limit  of  small  y.  We  then  get 

^  =  us(i(0,0  (29) 

~  =  Vu,(i(/),f)  y  (30) 

Equations  (29)  and  (30)  are  the  approximate  singularity  trajectory  equations  when  a  complex 
singularity  is  close  to  the  physical  domain.  In  this  case,  one  needs  xis  and  its  gradient  in  only 
the  real  domain,  which  can  be  determined  numerically. 

H’omplex  stable  manifold  of  a  real  fixed  point  of  (27)  is  defined  as  the  set  of  complex  initial  conditions 
/(,(())  such  that  satisfying  (27),  tends  to  the  fixed  point  as  l  — >  oo. 
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Once  the  singularity  locations  are  determined,  one  can  determine  p(:ro(f),  t)  and  Vd(:?o(/),  t ) 
in  terms  of  a  relatively  simple  system  of  ordinary  differential  equations,  as  will  be  derived 
shortly.  These  quantities  are  of  interest  since  for  a  point  x  sufficiently  close  to  To(0  (so  that 
d  is  small)  but  not  on  the  singularity  surface,  we  have 

u5  ~  a[Vd(x0(t),t)  ■  (x  -  x0(t))]-{'~a)p(x0{t),t)  (31) 

For  brevity  of  notation,  we  define 


and  (22)  implies  that 


A(t)  =  Vd(x  0(t),t) 

(32) 

B{t)  =  q(x0(t),t) 

(33) 

C(t)  =  p(x0(t),t) 

(34) 

t)  =  0  and  each  of  /,  f  /  f  are  zero  at  d  = 

0,  then  (21) 

AB  =  0 

(35) 

Ax  B  =  C 

(36) 

C  x  A  =  (A  ■  A)B 

(37) 

It  follows  that 


Furthermore,  it  is  clear  that  (25)  implies  that 

^  =  (us-A)B  +  TC 

where  T  is  a  second  rank  tensor  whose  elements  are  defined  by 

T,t  = 


Assuming  A  ■  A  /  0^  we  get  from  (37)  that 


f  - sc?  <40> 

where  S  is  a  second  rank  tensor  whose  elements  are  given  by 

— » 

Sjk  =  Tjk  +  tjkiAi  — y  (41) 

A  ■  A 

'Note  that  it  is  possible  for  A  ■  A  =  0  without  A  being  zero  since  the  dot  product,  defined  here  involves 
no  complex  conjugate.  However,  this  is  not  so  if  A  is  real  or  purely  imaginary. 
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given  that  tjki  is  the  usual  Levi-Civita  tensor  and  the  Einstein  summation  convention  on 
repeated  indices  has  been  used. 

We  now  argue  that  A  remains  nonzero  when  initially  nonzero.  By  taking  the  gradient  of 
(24),  it  is  clear  that 


dA 

dt 


-T +A 


(42) 


where  the  +  superscript  refers  to  the  transpose  of  T  defined  by  (39).  Therefore,  the  only 
way  for  A  to  be  zero  in  finite  time  is  if  the  elements  of  the  tensor  T  blow  up  in  finite  time. 
If  us(x,t)  is  smooth  for  all  times  and  grows  at  infinity  (for  an  unbounded  domain)  no  faster 
than  linearly  in  x ,  then  clearly  T(xo(t),t)  cannot  blow  up  in  finite  time.  Thus,  A  cannot 
be  zero  in  finite  time  when  initially  nonzero.  Equations  (40)  and  (42)  form  a  closed  system 
of  equations  for  A  and  (7,  and  thus  for  S7d(x0(t),t)  and  p(xo This  system  of  ordinary 
differential  equations  can  be  solved  numerically  once  us  and  its  spatial  derivatives  are  known 
in  the  complex  plane.  In  the  limiting  case  of  a  singularity  approaching  the  physical  domain, 
further  simplifications  are  possible  making  it  feasible  to  calculate  A  and  C  in  an  approximate 
manner  from  knowledge  of  us  and  its  gradient  in  the  real  domain  only. 

The  motion  of  complex  singularities  of  vorticity  is  significant  since  its  approach  towards 
the  physical  domain  can  be  associated  with  the  appearance  of  small  scale  structures  at 
places  in  the  real  domain  that  are  closest  to  the  complex  singularities.  Furthermore,  it  is 
clear  that  the  location  of  initial  singularities  in  the  complex  plane  is  related  to  the  data 
specification  in  the  real  domain  in  an  ill-posed  way;  i.e.,  one  can  specify  data  in  the  real 
domain  with  arbitrary  but  finite  precision,  yet  one  cannot  nail  down  the  singularity  location 
in  the  complex  spatial  domain  to  any  precision.  Yet,  if  the  complex  singularities  approach 
the  physical  domain  along  some  trajectories,  one  can  explain  the  sensitivity  of  the  later  time 
evolution  in  the  physical  domain  to  initial  conditions.  One  can  conceivably  start  with  a 
random  distribution  of  singularities  in  the  complex  domain  and  generate  an  understanding 
of  where  and  when  in  a  physical  flow  one  can  expect  fine  scale  structures  where  viscosity  is 
important.  The  full  time-dependent  calculation  will  be  the  subject  of  a  future  investigation. 


3.  POSSIBLE  IMPLICATIONS  FOR  THE  HYDRODYNAMIC  STABILITY 
OF  STEADY  FLOWS 


In  order  for  a  flow  to  be  stable,  any  disturbance  to  a  basic  steady  flow  cannot  grow 
without  bound  as  t  — *  oo.  A  flow  will  be  considered  absolutely  unstable  here  if  there  is  a 
continual  growth  of  disturbance  vorticity  at  some  fixed  Eulerian  point  in  the  flow  field.  It 
will  be  considered  unstable  in  a  Langrangian  sense  if  there  exists  a  disturbance  that  grows 
at  some  real  point  that  moves  with  the  steady  flow. 

If  we  consider  a  basic  steady  flow,  and  associate  the  steady  flow  velocity  and  vorticity  with 
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us(x)  and  u3s(.r),  respectively,  as  described  in  the  last  section,  we  can  consider  superposition 
of  an  initial  disturbance  so  that  the  initial  condition  is  of  the  form  (8)  and  (9).  Our  arguments 
in  the  last  section  suggest  that  any  point  fo(0  on  the  singularity  surface  d(x0(t),t )  =  0  will 
move  according  to  (27).  Since  us  is  now  a  steady  flow,  (27)  reduces  to 


dx  o 
dt 


=  us(x0(t)) 


(43) 


When  a  singularity  continually  moves  closer  towards  the  physical  domain,  one  can  expect 
smaller  and  smaller  scales  to  be  amplified  in  the  physical  domain  at  some  physical  point  in 
the  immediate  vicinity  of  the  complex  singularity.  More  precisely,  in  the  limit  of  a  complex 
singularity  approaching  the  real  domain,  its  distance  from  the  real  domain  will  be  of  the  same 
order  as  the  small  scales  generated  in  the  real  domain.  As  discussed  before,  singularities  of 
the  form  considered  here  can  at  best  approach  the  physical  domain  continually  without 
actually  impinging  it  in  finite  time. 

The  relevant  question,  as  far  as  hydrodynamic  stability  in  the  absolute  sense  is  concerned, 
is  whether  the  disturbance  vorticity  f'(d(x,t))p(x,t)  grows  indefinitely  with  t  at  some  real 
point  x.  Since  d  can  never  be  zero  in  the  real  domain  and  must  be  bounded  below  by  the 
smallest  value  of  d  on  the  real  domain  at  the  initial  time  (which  follows  from  (24)),  it  is  clear 
that  the  growth  of  disturbance  vorticity  will  depend  on  the  growth  of  p  in  the  real  domain. 

We  will  now  show  that  if  there  is  a  trajectory  £o(0  determined  from  (43)  that  asymptotes 
to  some  real  fixed  x s  as  t  —>  oo,  then  it  is  possible  to  choose  a  disturbance  such  that 
p(x0(t),t)  grows  exponentially  in  time. 

Without  any  loss  of  generality,  we  will  choose  xs  =  0.  Then,  according  to  the  given 
condition,  there  exists  an  x0(£)  such  that  it  approaches  the  origin  as  t  — ►  oo.  Near  the 
origin,  i.e.  the  stagnation  point,  it  is  clear  that 


us(x )  ~  Tx 


(44) 


where  T  is  now  a  constant  matrix.  Then  according  to  (43)  and  (44),  for  t  >  t0  (where  t0 
is  sufficiently  large),  x0(f)  is  determined  approximately  by 


(45) 


According  to  our  condition  of  approach  to  the  origin,  T  must  have  at  least  one  eigenvalue  a 
w.'n  a  negative  real  part.  We  choose  x0(t0)  to  be  in  the  eigenspace  corresponding  to  a.  If  a 
is  real,  x0(to)  will  be  chosen  to  be  completely  imaginary.  Whatever  the  case,  we  have 


x0  (t)  =  f0(*oK('~to) 


(46) 
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Then,  if  we  choose  A(to)  =  £  to  be  in  the  eigenspace  of  T+  corresponding  to  the  eigenvalue 
ex,  according  to  (42),  it  follows  that 

A(t)  =  £r<T,‘-<c)  (47) 

From  (46)  and  (47),  it  is  clear  that  A  ■  x0  is  bounded  as  £  — »  oo.  Equation  (40)  now  holds 
where  the  expression  (41)  for  Sjk  is  given  by 

Sjk  =  Tjk  +  fjt/6  -^-4  (48) 

Clearly,  since  £  is  a  left  eigenvector  of  T  (i.e.  an  eigenvector  of  T+),  it  is  cle<.  from  (48)  that 
it  is  also  a  left  eigenvector  of  S.  Thus  the  set  of  eigenvalues  of  T  .  nd  S  are  identical.  Since 
the  velocity  gradient  tensor  T  is  traceless  (incompressibility  constraint),  it  is  clear  that  the 
sum  of  eigenvalues  ol  T.  and  therefore  of  S,  must  be  0.  Thus,  since  there  is  an  eigenvalue  a 
with  a  strictly  negative  real  part,  there  must  exist  some  other  eigenvalue  A  with  a  positive 
real  part.  If  £  is  a  nonzero  vector  of  S  in  the  subspace  corresponding  to  eigenvalue  A,  it 
follows  that  if  C(£0)  —  £,  then 

C(t)  =  (49) 

Thus,  clearly  there  exist  disturbances  for  which  C  and  therefore  p{x0(t),t)  grows  exponen¬ 
tially  as  £  — >  oo.  Note  that  the  choice  made  in  (49)  is  perfectly  consistent  with  the  condition 
C ■  A  =  0  that  follows  from  (37),  since  £  and  £  are  left  and  right  eigenvectors  of  S  correspond¬ 
ing  to  differing  eigenvalues  a  and  A.  Our  proof  that  p(x0(t),t)  grows  exponentially  in  £  when 

— #  - * 

To (£)  approaches  xs  requires  that  £  •£  ^  0.  In  the  unusual  case  when  there  is  no  eigenvector 
of  T+  with  this  property,  a  modified  analysis  is  required  by  choosing  a  purely  imaginary 
vector  £  in  the  subspace  corresponding  to  the  eigenvalues  a  and  <r*  (Note  that  eigenvalues 
with  a  negative  real  part  must  be  complex  in  this  case).  Then,  the  tensor  elements  Sjk  in 
(48)  will  not  be  constants  but  are  periodic  in  time.  In  that  case,  we  have  been  unable  to 
show  in  any  elementary  way  that  C(t)  must  grow  with  time,  though  we  expect  this  to  be 
the  case. 

Note,  however,  that  for  two-dimensional  disturbances  superposed  on  '-dimensional 
flow  with  a  stagnation  point,  C(£o)  must  necessarily  be  aligned  along  the  x 3  axis,  which  is 
an  eigenvector  of  T  corresponding  to  eigenvalue  0;  therefore,  it  cannot  be  chosen  to  be  £. 
Indeed,  for  two  dimensional  disturbances,  the  terms  on  the  right  hand  side  of  (38),  which 
arise  from  vortex  stretching,  are  absent  and  C  is  a  constant. 

From  (24),  it  is  clear  that  at  the  stagnation  point  T,,  d(Ts,£)  =  d(xs,  0).  Thus  if  d(xs,  0) 
is  chosen  sufficiently  small,  each  of  /,  // /',  /2 f"j f*  appearing  in  (25)  is  going  to  be  uni¬ 
formly  small.  Thus,  it  seems  plausible  that  as  x0(t)  approaches  xs  and  /7(7o(£),£)  grows 
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exponentially  in  time,  p(xs,t)  (and  hence  the  disturbance  vorticity  f'(d)p(xsJ))  will  grow. 
However,  the  mathematical  proof  of  this  statement  is  far  from  obvious  for  a  general  flow, 
where  p  is  determined  by  the  nonlinear  equation  (25).  However,  in  the  Appendix,  under 
some  additional  assumptions  which  appear  reasonable,  we  are  able  to  show  that  exponential 
growth  of  p(£o(t),t )  does  correspond  to  a  growth  of  disturbance  vorticity  =  f'(d)p  at  xs 
for  a  steady  flow  where  Vms  is  a  constant  and  up  is  small  enough  to  permit  linearization. 
We  expect  that  the  connection  actually  transcends  these  restrictions,  though  we  have  not 
succeeded  in  proving  this. 

Based  on  the  arguments  in  this  section  and  the  results  in  the  Appendix,  we  present  the 
following  conjectures. 

Conjecture  1 

A  sufficient  condition  for  the  absolute  instability  of  a  smooth  steady  Euler  flow  in  a 
periodic  or  unbounded  domain  is  that  a  singularity  moving  in  the  complex  plane  according 
to  (27)  has  the  property  that  as  t  — >  oo.  xo(t)  — >  xs  for  a  real  xs ,  where  Im  Xo(0)  ^  0. 
This  condition  is  equivalent  to  the  existence  of  a  stagnation  point  xs,  where  the  velocity 
gradient  tensor  T  has  at  least  one  eigenvalue  with  a  negative  real  part.  In  that  case,  the 
vorticity  at  xs  will  grow  in  time. 

Remark  l 

The  above  conjecture  is  only  a  sufficient  condition  for  instability,  though  not  necessary, 
as  shown  in  a  later  example.  This  is  because  only  specific  forms  of  disturbances  have  been 
chosen-  namely,  those  associated  with  a  specific  singularity  structure  in  the  complex  plane. 
Remark  2 

There  is  no  restriction  on  the  size  of  the  disturbance,  though  there  is  a  restriction  on 
its  form.  Nonetheless,  we  do  not  expect  this  to  be  useful  in  uncovering  nonlinear  instability 
of  flows  which  are  linearly  stable  since  the  condition  for  approach  of  singularities  is  equally 
valid  for  linearized  equations. 

Remark  3 

Arbitrary  disturbances  in  the  conjecture  refer  to  any  arbitrary  disturbance  in  three  dimen¬ 
sions,  although  the  steady  flow  can  be  two-dimensional.  For  two-dimensional  disturbances 
superposed  on  two-dimensional  flow,  the  approach  of  a  complex  singularity  can  be  associ¬ 
ated  with  the  growth  of  the  gradient  of  vorticity  only.  This  can  be  seen  by  just  taking  the 
gradient  of  the  scalar  vorticity  equation  and  evaluating  it  at  the  stagnation  point,  which 
leads  to  dC/dt  =  TC,  where  C  —  Vt o(x3,t),  and  T  is  the  two-dimensional  velocity  gradient 
tensor  at  the  stagnation  point. 

Remark  4 

This  conjecture  can  be  applied  to  any  steady  state  flow  -  even  ones  that  are  numerically 
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determined  -  once  the  deformation  tensor  T  at  a  stagnation  point  is  obtained. 

Conjecture  2 

A  steady  flow  is  unstable  in  the  Langrangian  sense  if  there  exists  a  choice  of  initial 
conditions  i(0),  t/(0),  A(0)  and  (7(0)  such  that  the  solution  to  (29),  (30),  (40)  and  (42)  have 
the  feature  that  as  t  — *  oo,  y  — >  0^  and  C(t)  — >  oo. 

In  the  following  subsections,  we  present  examples  of  the  application  of  Conjecture  1  to 
specific  simple  Euler  flows.  Conjecture  2  appears  to  be  harder  to  apply;  at  this  point,  it  is 
not  clear  to  us  how  this  can  be  applied  in  a  fruitful  manner. 

3.1  Elliptical  Flows 

We  first  consider  the  case  of  plane  strain  in  an  infinite  flow  domain  with  the  velocity  field 


us(x,  0)  =  (ax\,  —ax2, 0) 


(50) 


with  a  >  0.  Then  according  to  (27),  a  singularity  at  Xq  =  (x\,x2,x3)  moves  according  to 


Hence 


dx  i 
dt 


OtX] 


dx2 
dt 
dx  3 
dt 


—ax2 


=  0 


x\  —  X|(0)ea<  ,  x2  =  X2(0)e  al,  x3  =  z3(0) 


(51) 

(52) 

(53) 

(54) 


from  which  it  is  clear  that  if  we  choose  Re  Xj( 0)  =  0  for  1  <  j  <  3  and  Im  xi(0)  = 
Im  ar3(0)  =  0,  then  x0(t)  — »  0  as  t  — >  oo.  According  to  Conjecture  1,  the  flow  is  going 
to  be  unstable  in  the  absolute  sense.  This  is  in  accordance  with  established  results  from 
hydrodynamic  stability  theory  (see  Pearson12,  Lagnado  et  al13  and  Bayley14). 

The  relation  of  the  approach  of  a  complex  singularity  to  the  appearance  of  increasingly 
small  scale  structures  is  highlighted  by  the  following  perturbation  of  our  straining  flow: 


u(x,  0)  =  (aa?!  -f  f(x2),  —ax2, 0)  (55) 

for  an  arbitrary  nontrivial  differentiable  function  /.  In  this  case,  there  is  a  simple  exact 
solution*’  given  by 

u(x,  t )  =  [axj  +  e~atf(x 2eat),  —ax 2, 0]  (56) 

^  Recall  *o(0  =  £(0  +  iy(() 

"This  is  probably  well  known;  however,  we  could  not  find  a  specific  reference 
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with  the  corresponding  two-dimensional  vorticity 


w  =  -f\x2eat) 


(57) 


Note  that  if  the  function  /  has  a  singularity  in  the  complex  plane,  it  approaches  the  physical 
domain  as  t  — >  oo.  For  instance  if 


fM  = 


(58) 


then  from  (55)  and  (56),  it  follows  that  the  singularity  of  u  and  u  is  located  at  x2  =  i  e~at 
which  approaches  zero  as  t  —*■  oo.  In  this  case,  with  a  two  dimensional  disturbance  super¬ 
posed  on  a  two  dimensional  flow,  the  approach  of  a  complex  singularity  does  not  correspond 
to  growth  of  vorticity;  however,  as  expected  from  Remark  5,  the  gradient  of  vorticity  grows 
exponentially  with  time.  The  solution  form  (57)  in  the  physical  domain  clearly  corresponds 
to  transfer  of  energy  to  smaller  scales  near  x2  =  0,  though  the  energy  contained  at  any 
specific  scale  eventually  decreases  with  time.  The  similarity  variable  eatx2  in  (56)  could  have 
been  predicted  by  looking  at  the  solution  (54)  in  the  complex  plane.  Note  that  Irn  Xj  and 
Iiii  x2  need  not  go  to  zero  for  all  initial  conditions  in  order  for  Conjecture  1  to  hold.  Indeed, 
there  exist  initial  conditions  when  singularities  move  further  away  from  the  real  domain 
which  corresponds  to  a  reverse  cascade:  energy  transfer  from  smaller  scales  to  larger  scales. 
This  can  also  be  seen  in  terms  of  exact  solutions  to  the  Euler  equation  of  the  form  (56)  when 
the  initial  perturbation  velocity  involves  the  second  component  rather  than  the  first. 

Three-dimensional  generalizations  of  the  above  straining  flow  can  be  made  for  which  the 
velocity  is  given  by 

us  =  (axi,  ftx2,ix3)  (59) 

where  a  +  /?  +  7  =  0  due  to  (11).  The  case  of  ft  =  7  =  —  corresponds  to  the  axisymmetric 
expansion/contraction  which  has  been  studied  extensively  in  the  turbulence  literature.  As 
7  — +  0,  the  plane  strain  case  considered  earlier  is  recovered.  Here,  we  show  how  to  apply 
Conjecture  1  in  its  simpler  form.  The  velocity  gradient  tensor  corresponding  to  (59)  is 


T  = 


/  a  0  0  \ 

0/9  0 
\  0  0  7  / 


(60) 


The  eigenvalues  of  (60)  are  a,  ft  and  7  -  at  least  one  of  which  is  negative  due  to  the  incom¬ 
pressibility  constraint  a  +  ft  +  7  =  0.  Hence,  according  to  Conjecture  I,  the  flow  is  unstable 
in  accordance  with  results  from  hydrodynamic  stability  (see  Pearson12  and  Townsend15).  We 
will  apply  Conjecture  1  to  the  remaining  elliptical  flows  to  be  considered  with  this  simple 
eigenvalue  test. 
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Plane  strain  plus  rotation  in  an  infinite  flow  domain  corresponds  to  the  base  flow 

us  =  (a^i  -  flxj  —  m'2,0)  (61) 

This  physically  corresponds  to  a  plane  strain  (with  strain  rate  a)  combined  with  a  solid  body 
rotation  (with  angular  velocity  ft)  in  the  plaine  of  the  strain.  The  velocity  gradient  tensor 
corresponding  to  this  flow  is 

/  cv  —ft  0  \ 

T  =  ft  -a  0  (62) 

0  0  0  / 

The  eigenvalues  of  (62)  are  0,  \/a'2  —  ft2  and  —  \/cv2  —  ft2  and,  hence,  one  of  them  has  a 
negative  real  part  if  cv2  >  ft2.  Consequently,  according  to  Conjecture  1,  this  flow  will  be 
absolutely  unstable  when 

a2 -ft3  >  0  (63) 

This  result  is  consistent  with  recent  stability  analysis  of  elliptical  flows  (see  Bayley14  and 
C'raik  and  Criminale16).  Since  Conjecture  1  only  provides  a  sufficient  condition  for  instability, 
it  does  not  formally  apply  when  ft2  >  cv2.  Nonetheless,  it  is  interesting  to  note  that  ft2  >  cv2 
is  actually  an  established  condition  for  the  stability  of  this  flow  (see  Bayley14). 

We  now  consider  the  base  flow 

us  —  (Sx2  —  20x2, 2ftaq,  0)  (64) 

which  serves  as  a  model  for  uniform  shear  flow  subject  to  a  system  rotation  with  angular 
velocity  ft  (see  Batchelor11).  Uniform  shear  flow  is  recovered  in  the  limit  ft  — ►  0.  The 
velocity  gradient  tensor  T  here  is  given  by 

/  0  S  -  2ft  0  \ 

T  =  2ft  0  0  (65) 

0  o  0 

The  eigenvalues  of  (65)  are  0,  ^2 ft ( .S'  —  2ft)  and  —  yj2$l(S  —  2 ft).  Hence,  according  to  Con¬ 
jecture  1,  this  flow  is  unstable  if 

2ft (.S’  -  2ft)  >  0  (66) 

since  this  condition  guarantess  the  existence  of  a  real  negative  eigenvalue.  In  terms  of  the 
rotational  Richardson  number  Ri  =  —2^  (l  —  2^  first  introduced  by  Bradshaw17,  it  is 
clear  that  according  to  Conjecture  1,  the  flow  is  unstable  if 

Ri  <  0  (67) 

This  result  is  identical  to  that  obtained  previously  by  more  complicated  linear  stability 
analyses  (see  Lezius18  and  Lezius  and  Johnston19). 
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For  a  uniform  shear  flow,  Ri  =  0  and,  therefore,  Conjecture  1  is  not  useful  in  determining 
the  stability  or  instability  of  the  flow.  The  same  is  true  for  the  more  general  cast-  of  plane 
shear  flows  where 

us  =  (f{x2),  0,0)  (68) 

Of  course,  it  is  well  known  that  uniform  shear  flow  is  weakly  unstable;  the  kinetic  energy  of 
infinitesimal  disturbance  grow  linearly  with  time13'15.  The  necessary  and  sufficient,  condition 
for  a  general  plane  shear  flow  (68)  to  be  unstable  is  for  /'(x2)  to  change  sign  -  the  classical 
inflexion  point  theorem  of  Rayleigh. 

3.2  Kelvin-Stuart  Cat’s  Eye 

This  class  of  exact  steady  solutions  to  the  Euler  equation  has  been  presented  by  Stuart20. 
Without  any  loss  of  generality,  his  solution  for  the  stream  function  can  be  written  as 


'|f(xJ,a:2)  =  -  In  [cosh  (2  x2)  +  p  cos  (2  Xi)] 


where  0  <  p  <  1  is  a  constant  characterizing  the  family  of  solutions.  The  velocity  field  us 
in  this  case  is  given  by 


sink  2x2 


p  sin  2xi 


\cosh  2x-2  +  p  cos  2xi  ’  cosh  2x2  -f  p  cos  2xi  ’  ) 

Within  a  period  7r  in  x\,  there  are  two  distinct  stagnation  points  at  (0,  0)  and  (|,  0).  At 
(|,0)  the  velocity  gradient  tensor  T  is  given  by 

(  0  A  °\ 


T=  \  -ri-,  o  o  (71) 

V  o  oo/ 

The  eigenvalues  of  this  tensor  are  0,  2ipx^2 /{\  —  p)  and  —2ipl^2/(l  —  p),  none  of  which  have 
a  negative  real  part.  Conjecture  1  is  therefore  not  applicable  and  nothing  can  be  deduced 
about  the  stability  from  the  existence  of  this  stagnation  point.  However,  at  the  stagnation 
point  (0,  0),  the  velocity  gradient  tensor  is 

(  °  k  »\ 

T  =  rfe  0  0  (72) 

0  0  0  / 


The  eigenvalues  of  this  tensor  are  0,  — 2p1/,2/(l  +  p)  and  2p1/,2/(l  T  p),  one  of  which  is  real 
and  negative.  Hence,  according  to  Conjecture  1,  the  Kelvin-Stuart  flat’s  Eye  is  unstable. 
This  agrees  with  the  findings  of  Pierrehumbert  and  Widnall21. 
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4.  CONCLUDING  REMARKS 

We  have  considered  the  motion  of  a  specific  type  of  singularities  of  the  incompressible 
three  dimensional  Euler  equation  in  the  complex  spatial  domain  that  correspond  to  a  smooth 
variation  in  the  real  domain.  Under  the  assumption  that  there  is  a  class  of  initial  conditions 
where  no  new  singularities  form  in  the  real  or  complex  domain  other  than  what  is  present 
initially,  we  have  derived  an  equation  (namely,  Eq.  (27))  for  the  trajectory  of  a  singularity 
and  have  shown  that  such  a  singularity  cannot  impinge  the  physical  domain  in  finite  time. 
Nonetheless,  we  expect  that  the  approach  of  these  singularities  towards  the  physical  domain 
will  correspond  to  the  formation  of  small  scale  structures  therein,  where  the  neglected  viscous 
effects  would  be  important.  The  sensitivity  of  these  small  scale  features  to  initial  conditions 
in  the  physical  domain  can  then  be  traced  to  the  problem  of  determination  of  the  initial 
singularity  location  in  the  complex  plane  from  finite  precision  data  in  the  real  domain. 
Furthermore,  from  (27),  we  note  that  the  problem  of  determining  the  singularity  trajectory 
is  the  complex  equivalent  of  determining  fluid  particle  trajectories  in  the  real  domain,  once 
an  Eulerian  flow  field  is  determined.  In  such  cases,  “Langrangian  chaos”  is  a  well-known 
phenomena.  Thus,  we  speculate  that  the  theory  of  dynamical  systems  may  be  relevant  to 
describing  chaotic  complex  singularity  trajectories,  which  could  have  a  bearing  on  turbulence. 

As  far  as  the  hydrodynamic  stability  of  steady  flows  is  concerned,  we  have  introduced 
two  new  conjectures  that  connect  the  complex  singularity  motion  with  the  growth  of  a 
disturbance  in  the  real  domain.  Admittedly,  more  theoretical  work  is  necessary  to  put  these 
conjectures  on  a  firm  foundation.  Nonetheless,  the  application  of  Conjecture  1  to  specific 
examples  that  include  the  class  of  elliptical  flows  and  the  Kelvin-Stuart  Cat’s  Eye  give 
results  in  agreement  with  previously  determined  results.  Though  no  new  instabilities  have 
been  discovered  so  far,  the  ease  and  simplicity  with  which  these  results  have  been  found  is 
certainly  encouraging.  Indeed,  as  mentioned  before,  Conjecture  1  can  even  be  be  applied  to 
a  numerically  determined  smooth  steady  Euler  flow  provided  there  is  at  least  one  stagnation 
point  where  the  velocity  gradient  tensor  has  an  eigenvalue  with  a  negative  real  part. 

So  far,  we  have  only  considered  the  stability  of  Euler  flows  with  no  boundaries.  We  do 
not  expect  boundaries  to  alter  the  statement  of  Conjecture  1  -  provided  that  the  stagnation 
point  is  in  the  interior  of  the  domain  -  since  the  instability  mechanism  associated  with  the 
approach  of  a  complex  singularity  is  associated  with  small  scale  disturbances  in  the  real 
domain.  This  issue  will  be  the  subject  of  a  future  study. 
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APPENDIX 


In  Section  3,  it  has  been  shown  that  if  the  velocity  gradient  tensor  T  at  a  stagnation 
point  xs  (taken  to  be  zero  without  any  loss  of  generality)  has  an  eigenvalue  a  with  a  negative 
real  part,  then  there  exists  a  singularity  trajectory  approaching  the  stagnation  point  such 
that,  for  sufficiently  large  f,  the  singularity  trajectory  x0(t)  asymptotes  to 

4(0  =  WMo)  MU 


In  that  case,  it  was  shown  that  p{x0(t0),  t0)  can  be  chosen  such  that  for  sufficiently  large  f, 
p(xo{t),t)  asymptotes  to 

p(x0(t)J.)  -  eHt~to)p(x0(to)Jo)  (A‘2) 

where  A  is  an  eigenvalue  of  T  with  positive  real  part. 

The  purpose  of  this  Appendix  is  to  show  that  the  exponential  growth  of  p  in  time  (as 
in  ( A 2 ) )  at  the  moving  singular  point  x0(t)  that  approaches  the  origin  (i.e.  the  stagnation 
point,  as  in  ( A 1 ) )  corresponds  to  a  growth  of  the  disturbance  vorticity  up  =  f'{d)p  at  the 
origin.  Here,  we  only  do  so  for  the  case  when  the  the  velocity  gradient  tensor  at  every  point  is 
T,  a  constant.  We  further  assume  that  the  disturbance  is  small  enough  so  that  the  following 
linearized  equations  hold 

+  u ,  •  Vt jd  +  uD  -  Vus  =  UD  ■  Vus  +  UJS  ■  WuD  (A3) 


V  x  U[)  =  uq 


where  =  f(d)  q. 

We  introduce  new  variables, 


V  •  ud  =  0 


xe-<T(t_(o)  =  y 


e°(l~to)CjD  =  W 


Then,  with  the  assumption  of  constancy  of  S7u,  =  T  (and  therefore  of  u;s),  we  have 
Wt  -  cW  -  aY  ■  VK W  +  TY  •  Vy  W  =  TW  +  u,  ■  VYuD 


Vy  X  Up  =  W 
V y  ■  =  0 

—4 

where  the  subscript  Y  denotes  a  derivative  with  respect  to  the  Y  components. 


Mio) 
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Equations  (A8)  -  (A  10)  are  linear  equations  with  no  explicit  time  dependence.  Assuming 
that  the  eigenvectors  corresponding  to  the  linear  stability  operator  form  a  complete  set, 
the  large  t  behavior  of  the  solution  for  a  generic  initial  condition  can  be  expected  to  be 
dominated  by  the  least  stable  (or  most  unstable)  eigenvector  times  cSt  (multiplied  at  best 
by  a  polynomial  in  t)  where  6  is  the  eigenvalue  with  the  largest  real  part.  The  analytical 
continuation  of  such  a  solution  to  the  complex  plane  will  have  similar  time  dependence.  Thus, 
according  to  this  argument,  the  time  dependence  of  IT’  at  a  fixed  complex  V  =  oxo(to)  can 
only  be  of  the  form  e6t  (multiplied  at  best  by  a  polynomial  in  /.). 

However,  for  the  specific  initial  condition  discussed  in  Section  3,  that  leads  to  (A2),  it 
follows  from  the  transformations  (A6)  and  (A7)  that  the  time  dependence  of  W(aio{to),t) 
is  of  the  form  efA+'7)<.  This  growth  rate  (or  decay  rate  if  Re  (A  +  a)  <  0))  cannot  exceed 
(or  have  a  decay  rate  less  than)  eSt  times  a  polynomial  in  t.  Thus,  Re  8  >  Re(  A  +  er) 
implying  that  for  a  generic  non-zero  choice  of  W( 0,  <o)i  the  growth  rate  of  IT'(0,  /0)  will  equal 
or  exceed  From  transformation  (A7),  it  follows  that  cZ?£>(0,f)  will  then  grow  at  a 

rate  at  least  equaling  eReSt,  implying  absolute  instability. 
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